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The ability to characterize static and time-dependent electric fields in situ is an important prereq¬ 
uisite for quantum-optics experiments with atoms close to surfaces. Especially in experiments which 
aim at coupling Rydberg atoms to the near field of superconducting circuits, the identification and 
subsequent elimination of sources of stray fields is crucial. We present a technique that allows the 
determination of stray-electric-field distributions (r), (r), (r)) at distances of less than 

2 mm from (cryogenic) surfaces using coherent Rydberg-Stark spectroscopy in a pulsed supersonic 
beam of metastable ls^2s^ ^Sq helium atoms. We demonstrate the capabilities of this technique 
by characterizing the electric stray field emanating from a structured superconducting surface. Ex¬ 
ploiting coherent population transfer with microwave radiation from a coplanar waveguide, the same 
technique allows the characterization of the microwave-field distribution above the surface. 


Hybrid systems aiming at coupling the internal state 
of atoms to solid-state devices have attracted significant 
interest in recent years. Realizations of such systems in¬ 
clude neutral atoms close to atom chips BHSl, tapered 
fibers 0 [TU] or photonic waveguides [TTJ |T2], Rydberg 
atoms close to mesoscopic devices [ME] and ions near 
surfaces [Mini. An important motivation for the de¬ 
velopment of hybrid systems is the combination of the 
long coherence times characteristic of atomic ensembles 
and the strong interactions and fast processing capabil¬ 
ities of solid-state devices as a route towards scalable 
quantum computing [IliMiEj- Atomic and solid-state 
systems can be coupled by electromagnetic fields which 
have to be controlled with high accuracy. For example, 
in our experiment we aim at realizing strong coupling 
between Rydberg atoms and solid-state circuit QED de¬ 
vices [M ED] using microwave fields. Stray fields ema¬ 
nating from the patterned surfaces of solid-state devices 
in general are of major concern in hybrid systems in¬ 
volving ions [TBl [T7] and Rydberg atoms [MIE| . The 
sources of stray fields are manifold and include surface 
adsorbates [MED], polycristalline surface patches [Mj - 
and charges in the isolating gaps of coplanar waveguides 
(CPW) [3S]. Although detrimental effects of electric 
stray fields might be mitigated in some cases, e.g., by 
microwave frequency dressing [55], coating the surfaces 
with adsorbates m or choosing chemically inert atoms 
(such as helium) [M [29], techniques to measure stray 
fields are essential. A possible technique in this context 
is based on Rydberg-electromagnetically-induced trans¬ 
parency, with which also the microwave field in a glass 
cell [30] or of a coplanar waveguide (CPW) [SU [32] has 
been characterized. 

In this article, we present a technique to measure static 
and time-dependent (microwave) electric fields above 
patterned surfaces in a cryogenic environment based on 
coherent Rydberg-Stark spectroscopy. We use a fast 
{v « 1700 m/s) supersonic beam of metastable ls^2s^ ^S'o 
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helium atoms. The atoms are excited to the 34s Rydberg 
state in two sequential, resonant one-photon transitions 
using a 10-ns-long laser pulse of wavelength A « 313 nm 
followed by a 160-ns-long microwave pulse of frequency 
zz « 27.966 GHz [55]. Specifically, we show that measur¬ 
ing spatially varying quadratic Stark shifts 

T ^ ( 1 ) 

= -Aa|i?(*)(r) -kE’^“'(rOP 

of the (conveniently chosen) 34s to 34p Rydberg-Rydberg 
transition for different well-defined applied electric fields 
[{Fx\r),Fj^\r),Fz''\f)), i = 1 — 4] allows the deter¬ 
mination of the three components of an unknown stray 
field (F®*''(r), emanating from the sur¬ 

face. Here, Aa = 1078.03 MHz (V/cm)“^ [MED] is the 
polarizability difference between the 34p and 34s states. 
The stray field is uniquely determined at each 

point r by a set of four independent Equations [Eq. Q , 
i = 1 — 4]. Three equations are needed because of the 
three dimensions of the problem and an additional equa¬ 
tion because of the quadratic nature of Eq. ([^. The 
spatially varying Stark shift is observed as 

the difference between the position-dependent resonance 
frequency r'a(r) of the 34s to 34p transition and its zero- 
field transition frequency vq = 27965.773 MHz, obtained 
from the known quantum defects [34]. 

We perform the measurements in an experimental re¬ 
gion cooled to ^ 4 K that is translationally invariant over 
6 mm along the helium-beam propagation direction (z di¬ 
rection) and bounded in the {x, y) plane by the surface of 
a patterned, 150-nm-thick superconducting NbTiN film 
on sapphire and a metallic shield (Fig.[^. The supercon¬ 
ducting film contains a coplanar waveguide with a 180- 
/rm-wide center conductor (aligned parallel to the z axis) 
separated from the ground plane by two 80-^m-wide insu¬ 
lating gaps. The chip is cleaned with an Ar'*' plasma prior 
to mounting under nitrogen atmosphere. Charges Qg and 
Qs (surface charge densities CTg and CTs, respectively) can 
accumulate in the insulating gaps of the waveguide [M 
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FIG. 1: Cross section of the interaction region above the 
coplanar waveguide. Potentials K, K, Vi and Vs can be ap¬ 
plied to the center conductor (green), the right (red) and the 
left (blue) ground plane of the coplanar waveguide (CPW) 
and to the metallic shield (brown). The thickness of the su¬ 
perconducting layer of the CPW is not to scale. Charges Qg 
and Qs (light and dark yellow) may accumulate either in the 
insulating gaps of the waveguide or on an insulating layer on 
the superconductor. The dashed, blue line indicates the cross 
section of the Rydberg beam limited by the collimating slit 
apertures. 


and on an insulating layer at the surface of the supercon¬ 
ductor. The insulating layer could originate from adsorp¬ 
tion of residual gases at the surface during cooldown [21] 
or from surface oxidation [35] . These charges represent a 
potential source of stray fields in our experimen¬ 

tal region. 

The region above the chip where the stray electric 
field is measured coincides in the {x,y) plane with the 
cross section of the atomic beam, limited by collimat¬ 
ing slits to the area marked by a dashed blue line in 
Fig.[T] and along x by the length of the Rydberg-atom 
cloud of 1 mm. Because of the translational invariance of 
the electrode configuration along z, Fz^\r) = 0, regard¬ 
less of the potentials applied to the device, see Fig. [^ 
Moreover, is constant and small enough that it 

does not cause shifts of the transition frequency of more 
than ~ 1 MHz, which is almost negligible compared to 
the shifts resulting from the x and y components of the 
stray field (see below). F^^'^{r) can be determined from 
the Stark shift at a position where the stray field in x 
and y directions has been compensated. Consequently, 
only the x and y components of the stray-electric-field 
distribution need to be determined, and only three sets 
of measurements {i = 1,2,3) in Eq. Q are required to 
determine F®*''(x,y) and F^*''^{x,y). 

Propagating along the CPW, the microwave square 
pulses of 200 ns length we use to monitor the 34s to 
34p transition have a Fourier-transform-limited band¬ 
width of ~ 4.4 MHz centered around a chosen detun¬ 
ing Aq from the field-free atomic resonance frequency vq. 
These pulses allow us to probe and image atoms subject 


to specific Stark shifts and thus specific field strengths, 
see Fig.j^a) for an example. This is achieved by ionizing 
the Rydberg atoms with a pulsed field of ~ 600 V/cm ap¬ 
plied 14 fj,s after the microwave pulse, when the Rydberg 
atoms have travelled more than 15 mm beyond the end of 
the chip, and extracting the electrons toward a particle 
detector [21]. The detector consists of two microchannel 
plates in chevron configuration and a phosphor screen, 
imaged by a CCD camera. An einzel lens operated at 
field strengths up to ^ 1.5 kV/cm expands the electron 
cloud by a factor more than 6 in x and y directions in 
a magnetic-field-free flight region of ~ 15 cm length, see 
suppl. material for more information. Because the 34p 
state lifetime is less than 2 ^s, only atoms in the 34s 
state are detected. The images therefore reveal regions 
of low intensities [dashed line in Fig. [^a)] with a spa¬ 
tial resolution varying between < 80 /rm and ~ 150 fj,ia 
which reflect the positions where atoms are efficiently 
transferred to the 34p state in the experimental region. 

Well-defined electric-field distributions are generated 
by applying potentials WiCpK- and 14 to the device 
(Fig.©. The applied potentials are small enough that 
they do not change the stray electric fields. The stray- 
field-measurement procedure outlined below starts once 
the stray fields have been roughly compensated by vary¬ 
ing the potential applied to the center conductor, a po¬ 
tential of 4 V being sufficient to reduce the Stark shifts to 
below 25 MHz over the entire experimental region marked 
by the dashed blue line in Fig. [T] 

We determine Ai5stark(2;, y)//i by varying the detun¬ 
ing Aq of the microwave pulses between —7 MHz and 
22 MHz and extracting the intensity at every pixel of 
the acquired camera images. This intensity corresponds 
to the relative population P^is{x,y) in the 34s Rydberg 
state after the microwave pulse. The amplitude of the 
microwave pulse is small enough to avoid coherent re¬ 
population of the 34s state. The spectra at each pixel 
are fitted by Gaussian functions to determine the width 
and detuning (i.e., the Stark shift) of the 34s to 34p 
transition. Fig. ib) shows, as examples, the spectra 
recorded for W = 4.00 V, V, = -0.02 V, = 0.00 V 
and 14 = —0.11 V at the two pixels marked by red and 
blue arrows in the image depicted in Fig. [^a). The fits 
yield the values AFstark//* = 0.15 MHz and 9.48 MHz, 
respectively, with corresponding widths of 4.2 MHz and 
14.67 MHz. The line broadening observed with increas¬ 
ing Stark shift is a consequence of the field inhomogeneity 
and the quadratic nature of the Stark effect at low fields, 
as explained in Fig. 4(a) of [55] . 

The measured distribution of Stark shifts [Fig. ic)] 
ranges from < 0.1 MHz (point marked by the label ’A’) 
to 22 MHz (below the point labeled ’+’), which is located 
close to the CPW. This range corresponds to electric-field 
strengths between < 16 mV/cm and 200 mV/cm, see 
Fig-i e). Because of the relationship between Stark shift 
and Stark broadening mentioned above, the observed dis¬ 
tribution of linewidths depicted in Fig. [^d) provides an 
independent measurement of the stray-field distribution. 
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FIG. 2: Distribution of electric-field strengths for potentials 14 = 4.00 V, V] = —0.02 V, 14 = 0.00 V and 14 = —0.11 V. (a) 
Measnred distribution of 34s Rydberg atoms after driving the 34s to 34p transition at a microwave detuning of 1.73 MHz. 
The spatial resolution of the images is 150 fim (see panels (c)-(f) for length scales). The dots mark the positions where the 
field-ionization signal is minimal and represent the positions where the 34s to 34p transfer is majcimal and hence the field 
strength (extracted from the Stark shift) is 55 mV/cm. (b) Experimental spectrum (dots) and Gaussian fit (solid line) of the 
34s to 34p transition for He atoms located at the positions marked by the blue (upper trace) and red (lower trace) arrows in 
(a) after binning over 5*5 pixels, (c) Spatial distribution of 34s to 34p transition frequencies (Stark shifts) extracted from 
the spectra measured at the different pixels, (d) Gorresponding distribution of full widths at half maximum, (e) Electric-field 
distribution derived from the observed Stark shifts, (f) Simulated electric-field distribution using Cg = —23.6(1) * 10“® G/m^ 
and (Ta = —2.10(5) * 10“® G/m^, see text. 


Over a large fraction of the cross section of the atomic 
beam, the observed microwave transitions exhibit a full- 
width-at-half-maximum of ~ 4.5 MHz [Fig. [^d)], close 
to the Fourier-transform limit of ~ 4.4 MHz of the mi¬ 
crowave pulse. This observation implies that the effects 
of inhomogeneous broadening that eventually limit coher¬ 
ence are negligible on the 200 ns timescale of the mea¬ 
surement. 

Using a 2 dimensional finite-element calculation, we 
simulate the electric-field distribution [Fig.j^f)] using the 
applied potentials 14, 14, 14 and 14 and fitting the sur¬ 
face charge densities Ug and (Tg (extracted best fit values 
given in the caption). The simulated field distribution is 
in good agreement with the experimental results and the 
comparison indicates that CTg and Us are defined with an 
accuracy of about 1%. Small discrepancies between the 
measured and simulated electric-field distributions result 
from the assumption that the charge density CTs is homo¬ 
geneous over the chip surface, which is the least strin¬ 
gent assumption we can make on the spatial distribution 
of these surface charges. The simulation also confirms 


the conversion factor of ~ 23 fj,m/px relating the camera 
pixels to the physical dimensions in the experimental re¬ 
gion, determined independently from simulations of the 
metastable He beam propagation and from simulations 
of the electron trajectories. From Fig. [^c-d), we extract 
an atom beam diameter of about 2.5 mm. 

To determine the components of the electric stray 
field Fy*''^{x,y)), we measure three different 

electric-field distributions {F^^f.{x,y), i = 1,2,3), for 
three different potential configurations [see Fig.[^a,b,c)]. 
The small homogeneous field component « 

50 mV/cm was determined from the Stark shift of ^ 
1 MHz at the field minima in Figs. [^a,b,c). 

(F®*''(a::,y),Fy*''(a;,y)) are extracted from a least- 
squares fit of the field strength |F(*^(a;,y) + F‘^^'^{x,y)\ 
to the measured electric-field strength F^^l{x,y) for i = 
1,2,3 using a random-search algorithm. The applied 
fields {x, y) were obtained from finite-element calcu¬ 
lations based on the applied voltages and the given ge¬ 
ometry. The reconstructed stray-field vectors displayed 
in Fig. I^d) all point toward the CPW, which reveals a 
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FIG. 3: (a-c) Distribution of stray-electric-field strengths extracted from Stark shifts measured for the potential configurations 
defined in each panel. The dotted line in (c) indicates a region of well-compensated stray fields (see text), (d) Electric-stray-field 
vectors extracted from measurements (a-c). (f) Predicted electric-field distribution as measured in panel (e), see text. Except 
for panel (d) all color scales are described by the color legend given in (a). 


negative potential difference between the CPW and the 
metallic shield and suggests an accumulation of negative 
charges (electrons) near the CPW. 

To test the validity of the reconstructed stray-field 
vectors [Fig. [^d)], a fourth potential configuration was 
explored, for which the field distribution depicted in 
Fig.i e) was measured. Comparison of this distribution 
with the field distribution displayed in Fig.J^f) obtained 
by adding the measured stray field [Fig. I^d)] and the 
calculated applied field indicates excellent agreement in 
the field strength, with deviations of at most 50 mV/cm 
in the region where the fields could be measured. 

The electric-field distribution depicted in Fig. [^c) cor¬ 
responds to a situation where the stray electric field 
is almost perfectly compensated over the broad region 
located above the red dotted line. In this region, 

\jFtot,Ax, y) + y) does not exceed 30 mV/cm 

after compensation, which corresponds to a reduction by 
a factor of more than 50 compared to the uncompensated 
stray-field distribution [Fig. [^d)[. 

A similar technique, based on measurements of the spa¬ 
tial variation of coherent population transfer, can be used 
to determine the distribution of the electric-field ampli¬ 
tude F^{f) of the microwave radiation generated by the 
CPW. Because of the translational invariance of the ex¬ 


perimental region along z, the rate with which atoms are 
transferred from the 34s to the 34p state (Rabi rate) only 
varies in the (x, ?/)-plane and is given by: 

\J ( 27 rA(a;, y)Y -b {dFt^{x, y)Y 
= \/ ( 27 rA(a;, y)f -b {s{x, y)Q)‘^. 

In Eq. 1^, A{x,y) = v — Vs,{x,y) is the detuning of the 
microwave frequency from the 34s to 34p resonance fre¬ 
quency Vs,{x,y), as determined from the stray-electric- 
field measurement, see Eq. 0. Because of the isotropy of 
the 34s state and the small static electric fields, the 34s to 
34p transition intensity (dipole moment d « 917eao) does 
not depend on the polarization of the microwave field. 
F^[x, y) is determined by measuring y) for a well- 

defined microwave amplitude 0 in the CPW and using 
the detuning A{x,y)^ known from the field distribution 
[see Fig. [^c)]. Using Eq. ([^, the microwave-field am¬ 
plitude is determined from the single parameter s{x,y), 
which scales the microwave amplitude in the CPW to its 
amplitude at the location (x,y). 

Fig. Qa) presents two measurements of the 34s to 34p 
population transfer carried out at the positions marked 
with arrows in Fig.j^b) by varying 0 from 0 V to 0max ~ 
560 p,V in the center conductor and using a 400-ns-long 
microwave pulse of frequency 27.966 GHz. Because the 
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stray-field distribution [Fig.j^c)], the displacement of the 
34s to 34p transition, and the lifetime of the 34p state are 
known, s{x,y) can be unambiguously determined using a 
simple model (black lines) based on Rabi’s formula (see 
Eq. (1) in The results are depicted in Fig. Bb)- 

The large dipole moment of the 34s to 34p transition 
implies that a microwave field amplitude as weak as ~ 
1.0 mV/cm suffices to fully transfer the population within 
400 ns. 

Overall, the measured microwave field decreases with 
increasing distance from the CPW, as expected. How¬ 
ever, the observed field maximum of the mode does 
not exactly lie above the CPW, but ~ 250 p,m to its 
left. Motivated by numerical simulations, this shift could 
be explained by interfering even and odd modes of the 
CPW because of different ground potentials of the (su- 
per)conducting surfaces with lateral dimensions larger 
than the microwave wavelength. 

In conclusion, we have presented a technique for mea¬ 
suring distributions of electric-field vectors with a spatial 
resolution of about 100 /im exploiting spatial variations 
of Rydberg Stark shifts and coherent population transfer 
in a supersonic beam of helium Rydberg atoms prop¬ 
agating above a cryogenic chip. By measuring several 
distributions of Stark shifts for the 34s to 34p transition 
resulting from different sets of potentials applied to the 
chip electrodes, the measured stray field could be accu¬ 
rately determined. The stray-field vectors allowed the 
identification of its source - in our case negative charges 
accumulating on the insulating surfaces in the vicinity of 
the CPW - and the determination of optimal potential 
configurations for stray-field compensation. 

Exploiting coherent population transfer between the 
34s and 34p Rydberg state, a similar technique allowed 
the measurement of the distribution of microwave ampli¬ 
tudes above the CPW. The reconstructed distribution of 
microwave field amplitudes did not exactly coincide with 
the expected mode of a CPW but revealed the effects 
of other interfering modes. By adding more microwave 
sources on the chip surface, it should be possible to fully 
characterize these spurious modes with the same tech¬ 
nique as used to determine the stray electric-vector field. 

For experimental setups that are not translationally 
invariant, it is straightforward to characterize the elec¬ 
tric field along the beam propagation axis (z direction). 
For different start times of the microwave pulse that are 
applied to the coplanar waveguide, the finite velocity of 
the atoms in z direction causes the atom cloud to be lo¬ 
calized at different positions along the beam propagation 
axis at which electric fields can be measured [29]. Also, 
the methods we presented are not limited to experiments 
with atomic beams but also applicable, e.g., to experi¬ 
ments using ensembles of ultracold atoms. 

The techniques presented in this article have several 
advantages compared to previous works using (Rydberg) 
atoms to determine electric fields. The data can be ac¬ 
quired rapidly, because only a single spectrum is enough 
to measure a full 2-dimensional electric field distribu- 
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FIG. 4: (a) Coherent population transfer between the 34s 
and 34p state as a function of CPW drive amplitude 0 for 
two spatially selected atom positions marked by the blue and 
red arrows in (b). The black lines are fits to a simple model, 
see text, (b) Microwave electric held distribution for ©max 
extracted from the measured distribution of Heff(rc, y) using 
Eq. §. The CPW is located at the origin (0, 0). 


tions, and a minimal set of 4 spectra characterizes all 
stray fields. Electric-field distributions are determined 
on a mm-sized area, limited only by the cross-section of 
the atomic beam. The spatial resolution of the Rydberg- 
atom images above the surface is currently limited to 
~ 100 fim in both transverse dimensions by the trans¬ 
verse velocity spread of the Rydberg atom beam. Simula¬ 
tions of the electron-trajectories through optimized elec¬ 
tron optics in combination with pulsed-field ionization 
indicate that laser cooling in the transverse dimension to 
the Doppler limit of ls^2s^ ^Si helium atoms (0.28 m/s) 
should enable an improvement of the resolution to 4 /rm. 
In contrast to other reported techniques, it does not al¬ 
ter the properties of the surface as it only requires low- 
frequency radiation (no laser radiation) close to the sur¬ 
face, and because of the use of an inert gas that does 
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not adsorb to the surface. Finally, it is compatible with 
cryogenic temperatures. 

Additionally, the spatial selection of atoms in field- 
free regions paves the way for precision measurements 
of atoms close to surfaces useful to study Van-der-Waals 
interactions m or surface ionization [38ll40) . Finally, 
the methods represent an important step toward cou¬ 
pling single Rydberg atoms to microwave photons on a 
chip. 
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